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ABSTRACT 


This thesis describes the develepment of a Fortran IV 
computer program for finding shock~induced stresses in a 
three-dimensional piping system. Discretized equations of 
motion are formed by the finite element method. input 
motions are support translations specified by response 
shock spectra. Maximum responses in individual modes and 
resulting octahedral shearing stresses at selected points 
are found for shock motions in three orthogonal directions. 
Latreme stresses, based on the square root of the sum of 
the squares of the modal stresses, are estimated for each 
input direction. An example problem is analyzed to demon- 


strate the use of the program. 
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described in Appendix E. 
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I. INTRODUCTION 


With the greater demand for electrical power, and the 
diminishing fossil fuel aipely: there has been an increase 
in the construction and use of nuclear power generating 
plants. Correspondingly, tnere has been greater concern for 
the safety of the power plant during ean earthauake. 

Alsc, the reduction in the size of the Navy has led to 
an increased emphesis on the integrity of the interna} 
systems o7 a naval vessel under shock loads. Thus, there 
has been a growing engineering interest in finding a rapid | 
and accurate method cf analysis to determine the shock- 
induced stresses in a complex continuous piping system. 

This thesis describes and presents 4 computer program, 
written in Fortran I computer language, to find the 
stresses in a piping system responding to shock. 

Complex piping systems can be analyzed by using dis- 
ecretized models. In using the discretizaticn techniaue, the 
continuous piping system is modeled by finite size elements 
connected at nodes. After the system has been subdivided, 
the elastic and inertial matrices of each element can be 
found and assembled to form the elastic and inertial 
matrices of the system [1,2]. 

A considerable volume of material is available to assist 
in analyzins structures and systems responding to 
earthquakes [3,4,5]. Two of the techniques used in earth- 


quake response analysis are modal analysis with a 


pine 


The 


discretized model and transfer functions [2,5]. ‘The shock 


input for eerthauekes can te svecified ay using tire 


history [3],or shock spectrum [4,51]. 


The Kavy uses @ techniau> celled the Dynemic Design-— 
Analysis Method (DDAM) to determine the shock response of 
shipboard equipment ene systems [6]. Tre DDAM is @ xodal 


analysis technigue using shcck srectr2 to specify the shock 
inputs. 

The shock spectrum vresents, @s a function of the systen 
natural frequency, the maximum response of a one Gegree-oFr-— 
freedom system responding to a shock motion. The shoek 
spectrum is generally ctresented in terms of resporse veiccity 
or acceleration. 

Several greduates of the Naval Postgraduate School have 
written theses on determining natural frequencies and mede 
shapes for piping systems. Fink [7] developed 2 prcgréenm to 
analyze planar piping systems, including cut-of-rlane 
bending, using transfer matrices. Baird [8] presented the 
necessary theory to expand Fink's work to a three-dimensicral 
piping system. Rudolf [9] developed a program that deter- 
mines the frequencies of a general three-dimensional piping 
system. because the transfer matrix does not yiel" a dis- 
cretized model of the structure, there is no ready means 
for using the frequency ard mode shape data to find shock- 


induceé modal responses. On the other hard, the finite 


element technique provides a consistent discretized medel 
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Here, the double errox«nheaaG vectors are positive rotations 


Tt 
a = [u, Uy --- YW ol 
From Fig. 1, components u, and uz are longitudinal dis- 


matrix with Frzemieniecki's transverse sheer deflection 


; i (- ee : 
Sion for the element xinetic energy 7” is formed, rotatory 


Sf ge gs lt 5 ae ae 
me = = Cz SQ UN NG dé = 24 oh G 
2 = eS et = = 


where £ is the element length, C is a constant, oe is the 


elemen® inertia ratrix, N is the shape function vow vector, 


= Ste E “ as poe f i : SS : ey. 
Therefore 
i 
e i ea Z 
Mo = C2 f N N dE (1) 
0 : 
Since the local element axes coincide with the cross 


section principei axes, the inertia matrix may be deter- 


and C = pA where A is the cross-sectional area of the 


element and p is the density. From Ea. 1 


b) twist: 


= Pow 8 bs 
Gm “ Cu, U6 9 Ny is (1 = E)g£, Nio = EL 


—_— 


and C = J where J is the mass moment cf inertia per unit 


length atcut the pipe «xis. From Eq. 1 


ue k 4 | (1b) 


c) bending in the 1-2 plane: 


: . Ah eae 
Gio = [uy Ug Ug Uy), Ny = 1-3e° ter, 


Ne = Ce-2e%4e3), Np = 3e°-209, 1, = a(-e%4e3), 


and C = pA. From Ea. 1 
156 
e _ phe 
Noi2 = woo | 24% 
5h 
-132 


ient se ; 
cients as B12 
rotation angles are reversed, 
and 4 are multiplied by minus 
assembled to form the element 


is displayed in Appendix A. 


. a STN = - <a 
279 54 1324 
he? 432 ~32° (1c) 
132 156 -229 
2 2 
—32° ~-222 He 


in the same numerical coeffi- 
positive senses for the 

rows 2 and # and columns 2 
one. These submatrices are 


inertia matrix. This matrix 


C. ASSEMBLED EQUATIONS OF MOTION 


The governing equations of motion, in matrix notation, 


for undamped free vibration 


[2] with a fixed base are 


(2) 


where q is the assembled displacement vector and MH and K 


are the assembled inertial and stiffness matrices. 


MH and 


K are referred to a global coordinate system and we and Ke 


are referred to local ccordinates. 


The assembly process 


consists of realigning the local system to the global 


system and building the ™ and 


K matrices. Reference [2] 


gives sample assembly techniques. 


te 


D. MODAL ANALYSIS 

Modal analysis is the process of finding the response 
motion of 2a many degree-of-freedom system by superposition 
of the response motions in the principal (or natural) modes. 
This techniaue is advantageous because the motions in the 
principal modes are independent of one another (uncoupled). 
To use the method, it is first necessary to find the nat- 
ural frequencies and mode shapes for the free motions 
governed by Ea. 2, i.e., to solve the eigenvalue problem 
Kv =a° Hy. 

Finding eigenvalues (w°) and eigenvectors (v) is a 
standard problem of numerical aralysis. The first ctep is 
a transformation to a single-matrix problem. Details of 
the transformation used here, which requires a Cholesky 
decomposition of the inertia matrix, are given on p. 349 
of reference [1]. Following this, Jacobi rotations [15? 
are used to find the spectral and modal matrices (2° and V). 
The modal matrix is normalized with respect to the inertia 


T 


matrix, i.e., Vo MV = 1, where £ is the nth order identity 


matrix. 

The shock analysis is based on using velocity shock 
spectra to characterize the input motion and mode parti- 
cipaticn factors to determine the responses. Although both 


shock spectra and mode participaticn factors have been 


extensively used in earlier studies [6], the finite element 


19 


discretization necessitates a new development of working 
formulas. Details of this development are given in 
Appendix B. 

Each shock input motion consists of base translation 
parallel to a global axis. It is shown in Appendix B that 


the mode participation factor by for mode vr is 


bo =v) wu (B.5) 


(r) 


where v is the modal eigenvector and u is the vector of 
absolute displacements corresponding to unit base displace- 
ment. The maximum relative displacements due to mede r 


response are 
sigh?) ay 
Vv Di Vis Os, (B. 6a) 


where ws is the modal circular frequency and Wi is the 
spectrum velocity at that frequency. 

Since responses to shock inputs in each of the three 
global directions are separately determined, the vector u 
and the mode participation factors by must be evaluated 
separately for each direction. The modal masses T,, (see 


Appendix B) are likewise different for each input direction. 


E. GENERALIZED FORCES AND STRESSES 
The subvector a° of element peak modal displacements 
may be found from the system vector, taking into account the 


eee 
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orientation of the local reference axes. From this the 


element generalized force vector ae is determined. This 


force vectcr consists of the elastic 


contributions tor each element. 


e (x® ae e) q® 


Pig. 2 shows the positive directions 


of the generalized force vector. 


Ge 
f — [f, f, eee 


and inertial 


(3) 


for the components 


FIG. 2 


ELEMENT GENERALIZED FORCES 
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The stresses at the nodal points of the element can be 
determined from the generalized forces by standard methods 


of solid mechanics [10]. 
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TIT. PROGRAM DEVELOPHENT 


The computer program was written in Fortran IY com 
puter language using double-precision arithmetic on an 
IBM 360/67 digital computer. It consists of a main calling 
program and seven subroutines. The main program calls the 
subroutines, where the actual calculations occur, in the 
necessery order. Information is passed between the main 
program and the subroutines by the use of a common storage 


core. Appendix E contains the program listing. 


A. SYSTEM DESCRIPTION 

ifter the piping system has been modeled and subdi- 
vided, subroutine INPUT is used to read the data input of 
the geometry and material properties of the system. The 
geometry of the piping system model is determined by the 
globeéi coordinates or the nedes of the model. The global 
system is a right-handed orthogenal coordinate system. 
Each pipe segment must be parallel to a global axis, but 
there are no restrictions on the placement of the origin of 
the global system. The parameters from the subdivision of 
the system are read. For each node the global coordinates 
and boundary conditions are read. If the node is fixed, 
the stiffness and inertial matrix components due to the 


node displacements are not assembled into the system 


matrices. 


Tne numbering cf nodes ane elements is arbitrary. 


Tnere are two nodes per element and six degrees of freedom 
per node. For each element, the ncde numbers and pipe 
group ere read. A pipe group includes all elements thet 
have the same Young's modulus, Poisson's ratio, specific 
weight, radius of gyration, and cross-sectional dimensions. 
Appendix D shows vJhe method used to calculate the effect 
of lagging on the specific weight and the radius of 


gyration. 


B. STIFFNESS AND INERTIAL MATRIX FORMULATION 

Subroutine F@RM generates the element stiffness ané 
inertia matrices as shown in Appendix A. Rotatory inertia 
and shear deformation due to bending are neglected. The 
orientation of the element is determined with respect to 
the global coordinate system and tne correspondence between 
local and global degrees of freedom is established. After 
determining the correspondence between the local and global 
degrees of freedom, the element stiffness and inertia 
matrices ere assembled to form the system stiffness and 
inertia matrices. If pipe hangers are present, tne corres- 
ponding (uniaxial) hanger stiffnesses are added to the 
appropriate diagonal elements of the system stiffness 


matrix. 


C. FREE VIBRATION MODE SHAPES AND FREQUENCIES 
Subroutine CiHAMOD uses Cholesky decomposition of the 


system inertia matrix and coordinate transformation to form 
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a singie syrcetric metrix for which eigenvalues amd cigen- 


vectors are Pound. Subreutine JACREY is a modification of ; 


@ program criginally coded by Professor @. Cantin. Tt 


Cowen Sd 


= 2 


uses the Jacobi variable threshoid methad te find eisgen- 


values and eigenvectors of @ re@i symmetric matrix. The 

model matrix, when transformed beck to system coordinstes, 

is. normelizeé with respect to the system inertie rmetriz. 
Sinee the systern natrices inere2se in size with Piner 


subdivisions cr more cco=vlex systers, En eronomtizins tech-— 

nigue, to reduce the number oF Gegrees af Freedom in the 

eigenvalue problem, was investigated Guring the progrem 

deveicrment. The component mece synthesis methed as used 

by Benfield and Zruda [11] wes stuGied. it was establisheé 

that the constrained component brench teechniave £113, «ith i 
interface loading, with no ncde suppression gaye resuits 

identical with those obtainec by the direct method descrired 

apove when anplied to planar vibration of 2 clamrped-cianped 

beam represented by six elements. Despite the success of 
this trial, it was conciuded that the added progrem ccm- 
plexity accomranying tne use of component mode synthesis 
vould offset the potential econcmies. Another standera 


economizer technique [14] was considered, but ultimately 


rejected for the same reason. 


D. MODAL RESPONSE 


The number of modes to te us- i in the stress analysis 


and the specification of the shocx spectrum are input data 


-“ 
ra 


for subroutine SPCTRE. FYhere 2re tire cholces arvatilanlie 


to select the mmber of medées to te used: 2 designated g 
tunner freauemcy limit, em esper Limit weich is a mmItiple 

of the fumcenemtal Freauentey, ami am ontion For any method 

the user desires. Three choices ere aisc exrailebie far 

specifying the stork suectrim: a smectrum defined by 

Straight-line segmerts; a2 comstemt velocity, them comsteant 
a@cesieraticn shock spectrum; end em ention for any method 


the user Gesires. Tre snrock imowt is restricted to a 
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Subrovtine STRESS caiculates stresses at th 
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the elements. The element peak Gisplacezents are ceter- 
minéa and realigned from giodal degrees of freedor beck to 
the local element degrees cf freedom. The element gener- 
alized force vector is then calculetecd. Censider the 


action of the generalizec force vector on the ¢lenent 


shown in Fig. 3. 


STRESS LOCATICN PORTS 


At icecations 2,5,c, and Gd, the stress vector ecting on 
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Tmese moma and simerirs stressss are them esed to Fim 


the octet=sdral sheerine stresses [121], weieh are given 
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printes fer the three shccx input directions. 
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iV. FROGRSN Testes 


Em order to explere the cencoilities of the sreogrez 
ama weri2y fits integrity, a mumfer of pleme test cnfiger— 
ations were stm@fed. Timse are ceserffed in the folloving 
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progr::.. The cresent rregrem gave identical resu 


mae ows 


39 


im-plans mo¢es. The cot-of-plame modes exhibited the 


reamired symmetric o> enti-spmmetric form. 

&gastions?! merogrem tests om this comfigmretion involved 
aifferest element mmmberins amd mode numberinms ami six 
aivverent oriemteations of the slahel ares retative to the 
structure. Mote shares ami freanemcies were mmefPected by 


ehese changes. 


B. YWES¥ PROSLEE 2 
This system, show: im Fic. 6, comsisted of 2 -treickt 


uniform pine clemped at the emis amd renresented by four 


elerents. 
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r 
solution with exact results [23]. All the eigenvectors 
showed the requirec symmetric or enti-symmetric forn. 
Bending moces occurred in péirs having equal eigenvalues 
and with the correspording eigenvectors representing 


deflections in orthogonal planes. 
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Stresses found for this system were studied in detail. 


ata 


She we nn 


At each mode, the stresses are caleulated at tuo points of 


the emoss—-section. 


2 


at these sectioms are caileulzted twice. 


wes found Detvecen ¢ 


Except at the clamped ends, there are 


two adjacent elememts sharing ceck mode, so thet stresses 


Comiete agreement 


mese XO SETS OF stress eveluetions. 


Also, identical stresses resulted from shocks im the tuo 


Girections perpemtiicular to the lensth. 


COMPARISON OF FERSY TERRES 
PaPhG S¥ST=M 2 WITH EXECT 


Mode Exact 
(22) 

i 822.32 
2 1364.35 
3 2262.20 


C. TEST PROBLEX 3 


Finite 

Element 
(Hz) 
£22.65 

4374.26 


2320.68 


SENDING 


FREQUENCIES OF TEST 
RESULTS 


Percent 
DifFerence 


the equal-legged ccenfiguraticn of Fig. 7 consists of 


uniform pive throushout. 


Mode shapes of this system egain 


showed the required symmetric or anti-symmetric forn. 


Likewise, the shock-induced stresses exhibited the expected 


symmetry properties. 
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V. EXAMPLE PROBLEM 


A. PIPING SYSTEM 

The exemple system analyzed is a model of the main 
steam piping system on a modern naval vessel. The main 
steam piping from the boiler 0 the rigid cross-connect 
anchor assembly is modeled. The ends of the piping system 
are clamped, and there are two hangers. The system is 


three-dimensional with nc branches. The pipe properties 


are: 
outside diameter 7.625 inches 
inside diameter 6.011 inches 
Poisson's ratio 0.300 
Young's modulus 30x10° psi 
specific weight 0.327 1b/in3 


radius of gyration 4,04 inches 


The specific weight is a fictitious value which accounts 


for the weight of five inches of lagging bonded to the pipe. 


Fig. € is a schematic of the piping system. The piping 
section Pl is 84 inches long, P2 is 264 inches long, P3 is 


108 inches long, and P4 is 288 inches long. 


B. DATA INPUT AND OUTPUT 

The input data cards are punched in accordance with the 
instructions contained in Appendix E. The echo-check of 
the data is shown in Table II and is arranged as it would 


appear at the end of the program deck. The system 
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FIG. 8 
EXAMPLE PIPING SYSTEM 


eigenvalues and eigenvectors are listed in Appendix C, 
Table VIIZ., Tables III and IV show the modal velocities 
and the mode particination factors. The stress outout is 
given in Table V. 

This problem required a storage capacity of 174K bytes, 


and 2.67 minutes of computer time. 
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MODE PARTICIPATION FACTORS, EXAMPLE PROBLE? 
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. FABLE Vo OCTASDREL SE=LAINS GYRESSES EreeTES FROSLEN, ESE : 


im -—-—— ; 


ELEMERT i-> + ELEMENT T : 


: 9.271250 03 6.523938 GS @.205355 O35 2.055150 GS 3.15562) 83 7.845550 63 t 
62674260 03 12687520 04 £.482565 TZ 42293349 03 5.395323 JS 2.8923£% C35 
26959595 G3 6.198625 33 72-2427210 O23 S.S26278 C3 3.455850 02 "32809223 63 

: 40588579 O35 82645532 03 3.995535 SS 42372620 63 5.55352 63 3.85000 03 * 


e 


ELEMENT = 2 
6.605590 03 6.392529 03 2.617210 
* 4588570 03 6.655522 03 3.095510 
1.879580 04 72527779 23 7.714059 0 
42960030 03 7392369 03 1.943069 0 


ELEMENT 6B : 
5.025879 33 32149459 23 3.869283 03 
4.372629 03 52550529 03 3.528030 63 
2.655362 03 3.213459 03 3.567759 23 
7242439) 03 1.318909 35 4.405299 33 


a oe mo 


WwW 


wv 


ELEMENT 3 7 LEHENT 9 
5.764590 02 72612435 C3 1.776203 03 72592123 03 24272929 0% 42365372 05 
1.888130 04 4.471280 93 7-52339D 903 22469200 03 2.611243 03 3.5989%D 03 
4.45858D 93 §.76172D 93 3.033889 O35 52776449 332° 1.013420 04 3.202215 02 
1.11412D 04 4.91465) 03 2.673222 2 857275 03 2.655160 032 8.55153 05 


ti 


. ZLEMENT = & , ELEMENT 10 
4.45358N 03 5.761720 33 3.03388D 93 54776440 03 1.019420 04 3.261110 03 
1611412D 04 4.91462D 03 2.673220 63 4.657270 03 2.063109 63 6.551305 03 - 
4.36188D 03 9.148555 03 3.568615 03 7.218965 G3- 6128639 03 4.520539 02 
754994BD 03 4.30096D 93 $.18924D 92 9.86817) 63 3.488539 03 1.570740 6% 
‘ELEMENT = 5 ELEHENT 11 
4.561580 03 9.146560 03 3.568610 03 6.540949 03 6.053360 93 4.586799 03 
7.499480 03 4.300965 63 6.185240 03 9-86817D 93 3.488530 93 1.579740 a4 
4454386D 03 7.59655) 33 4%.76524D 03 1.043000 04 4.379220 03 2.33575D 63 
6.252800 03 4.33059D 03 7.981600 03 3453535D 03 2.30537D 03 9.55294D 03 


ELEHENT 6 - ELEAENT i2 


4.54386D 03 7.596560 032 4.765240 02 1.943930 04 4.379220 9) 2.335750 93 


2 
6.252800 03 %.33659D 93 7.981600 93 3.535340 03 2.395370 02 9.55394D 03 
136613900 03 4.321470 03) 2.472250 03 1.769770 04 1.275180 04 6.00477D 03 
1,C386730 04 4.381410 03 7.711135 03 6.153500 03 2.751590 33 1.630620 04 


VE. DESCUSSIGN 
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Some of tne Linftetioms of the cenebilit2. 


progrem develoned above are comsidered here. 


&. CONFIGURATION LMHPraTEONsS 


can be modeled. Fhere are no imrerent iimitetionse om aliexw- 
abie tonologic si corplexity. 


are not represented in the present modeline include edded 
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For the present purpcse, the appropriate measure of 
system size is the numher of degrees of freedom of the 
system model. This number n, which bears no direct relation 
to physical size, determines the core storage requirements 
and execution time of the program. For approximate estima- 


2 
tion, storage requirements are proportional to n° and 


3 


execution time is proporticnal to n-. Using the date from 


tne example problem, one can estimate that a 100 
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éeorec—of-—frecdam system world recuire shout OTE bytes oF 
eGorte Statage ent ehact 9 mimntes execution time. 

Trerecesing the orablem sire aise imerceses tim mommad-—ofr 
errezs. Becemse dacdle-orecisian eritimetic (55 bit memtissa} 
is msed throusmont, atserved roumti-off effects have teen 
founm® meslicthle (2 meaztmem of 7 memfit im the Sth signifficant 
Gisit Zor m = £5) im aki anolicetions to dete. The very 


smal? errams Getected are attributed to the residual eicen- 


terminated. 
Hm viexu oF the foregoins comsiderations, it is believed 


eat, Fa Saks pea as e Rees 
thet the prectical upper limit on protiem size (with the 


3BH 360/67) is about m = 130 en€ is determined by core 


storage cavecity. 
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He is conmcinfed that am effective gresrem hes teen 
@eveioped for cstermining shock-indreed stresses in pipims 
svstems. The orimeinal Minitetions oF the present versio: 
een be remoret Dy etdims Features that are clearly within 


ys ~ - % ~— ne ee 7 ~ 
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~, =* = = to = a 
i. Remove the gine exes grientetion restricticn so that 
= z — —— ~, ni 
@ gencrei pinpims Systes cam te analyzed. 


—J a = Om, ~ RJ = oo os - 
2. Develico slement stiffness emi inertia metrices for 


= + rt o 4 a—sttraqz Ry . L235 as = 
5. Repiece the dJacoti retetion mevhoc Sy 2 more efricient 
eigenyalue-cizenrector 2igerithm. 
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APPENDEX B=: 


VWEROQCETYS SHOCK SFECTRA 20D MODEL RESPONSS” 


A shock srectrum ezmibits the maxim resvomse dis-— 


plzcement of a simsie degrec—of-—freedom system whose bese 


32s subjecced to the shock motion. 


Consider a sinsie degree— 


s-freetom whose these hes 2 shock displacement Ss, and uhose 


Gisplacesent reletive to tm bese is Z Tne 


eanetion of 


treme value of z in response tc the shock motion s, then the 


Gisplecerment sncck spectrum for that shock motion is 2 plot 


(B.2) 


To develop the eauations for finding modal response 


from shock spectrum data, let w be the vector of absolute 


1 


This Appendix is based on material presented ty Prof. 


R.E. Newton in tne covrse NE4522, September 1971 . 
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Gisplecements im the system desrees-—of-freedom. For 2 


umiaxiai tramslation s of the sese, this my be exnressed 
as 
w= @tus (-3) 
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where gis the vector of é 
base, end u is a veetor of constants. To illustrate the 3 


meanings of these vectors, comsider the bees element of 


Acie. Gael a 


Base 


BEAM ELEMENT COORDINATES : 3 


: = T 7 T 
For this example, q = [a, 4, a, Q,)]°, u = [a,+#s a, a,*s dy], 


and u = {10610]. It can be seen that u represents the 


ae a eta 6 
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thee 
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v x 


‘me equations of motion of an nm degree—of-—freedon 


system with base motion may be written 


Muwet+Ka=0 


Using the substitution g = V p where p is the vector of 


principal ccordinetes, this may be rewritten as 


Pp oo rt m od 
SN a Ree ane 
or 
- 2 . oo 
where b = vi Mu. Eq. B.4 is equivalent to the scalar 
equations 
D,+tu,?p,=-b, 8, (re1,2,...,n) (B.4a) 
r r Pr cr > rane 8 ad ‘ 


where p,, and b,, are the rth components of p and b, 


respectively, and w, is the circular frenuency cf mode r. 
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3 
fhe component bd, is called the mode participation factor ; 
for mode r. It is given by i 

b= vw (B.5) 


where vy) is the eigenvector for mode r {rth column of 
v). 
Comparing Eq. B.4a with Ea. B.1, and using Ea. B.2, 


the maximum response in mode r may be expressed as 


=hb v Sos £ 
(D.) max by V)/e, (B.6) 
where V, is the spectrum velocity at frequency w,° Using 
Eq. B.6, the corresponding extremam for the relative 


displacement vector q is 
ev) yo Vy (B. 6a) 
Gmax ~ © py Si 


Eq. B.6a is used for calculating peak modal responses in 
order to find stresses. 

In shipboard shock studies, it has been found that 
massive equipment may partially suppress the base mction 
[6]. This may be taken into account by appropriate modifi- 
cation of the spectrum velocity Vis Such a correction re- 
quires apportionment of the system mass among the modes. 


The modal mass for the rth mode is taken to be 


fo 2 
m, = DL (B.7) 
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APPENDIX C: EXAMPLE PROBLEM EIGENVALUES AND EIGENVECTORS 


Table VIII shows the first six eigenvalues and the 
corresponding eigenvectors of the example problem. The 
eigenvalue (square of the circular frequency) appears at 
the head of a column followed by the components of the 
eigenvector grouped by nodes. The remaining eigenvalues 
and eigenvectors are omitted. Spacing of the eigenvalues 


remains approximately uniform throughout the remainder of 


the set. The largest eigenvalue is 1.32 x 10°, 
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APPENDIX D: CALCULATION OF MODIFIED SPECIFIC WEIGHT j 


AND RADIUS OF GYRATION 


There is negligibie effect on the stiffness of the 
pipe element from the lagging on the pipe; however, the 
lagging dses contribute significantly to the mass of the 
element. By appropriately modifying the pipe element 
specific weight and radius of gyration, the mass effect of 
the lagging can be included. The modified specific weight 


is 
Yq = (Wy + Wy )/A (C.1) 


where en is the modificd specific weight of the pipe ele- 
ment, and the sum W, + Wy) is the combined weight of the 
lagging and pipe per unit length. The radius of gyration 
is given by the relation 


% 2 ; 
a + J. = (m, + m,) r (C.2) 


L 
where Jy and Jy, are the mass moment of inertia per unit 

length of the pipe and lagging, respectively, (m, + ny) is 
the combined mass of the pipe and lagging per unit length, 


and r is the radius of gyration. 
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APPENDIX E: PROGRAM LISTING 
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